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Abstract 

This paper provides an introduction to some stochastic models of lattice gases out 
of equilibrium and a discussion of results of various kinds obtained in recent years. 
Although these models are different in their microscopic features, a unified picture 
is emerging at the macroscopic level, applicable, in our view, to real phenomena 
where diffusion is the dominating physical mechanism. We rely mainly on an ap- 
proach developed by the authors based on the study of dynamical large fluctuations 
in stationary states of open systems. The outcome of this approach is a theory 
connecting the non equilibrium thermodynamics to the transport coefficients via a 
variational principle. This leads ultimately to a functional derivative equation of 
Hamilton-Jacobi type for the non equilibrium free energy in which local thermo- 
dynamic variables are the independent arguments. In the first part of the paper 
we give a detailed introduction to the microscopic dynamics considered, while the 
second part, devoted to the macroscopic properties, illustrates many consequences 
of the Hamilton-Jacobi equation. In both parts several novelties are included. 
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1 Introduction 



Models have played a fundamental role in equilibrium statistical mechanics. The Ising 
model provided the first proof that statistical mechanics can explain the existence of 
phase transitions and was a main guide in the study of critical behavior. A reason for this 
effectiveness is the circumstance that the macroscopic behavior is, to a considerable extent, 
independent of the microscopic details. Hence different systems exhibit qualitatively the 
same phenomenology at large scales. 

Out of equilibrium the situation is more complex. First the variety of non equilibrium 
phenomena one can conceive makes it more difficult to define general classes of phenomena 
for which a unified study is possible. Furthermore the details of the microscopic dynamics 
play a far greater role than in equilibrium. Since the first attempts to construct a non 
equilibrium thermodynamics, a guiding idea has been that of local equilibrium. This 
means the following. Locally on the macroscopic scale it is possible to define thermody- 
namic variables like density, temperature, chemical potentials... which vary smoothly on 
the same scale. Microscopically this implies that the system reaches local equilibrium in a 
time which is short compared to the times typical of macroscopic evolutions as described 
for example by hydrodynamic equations. So what characterizes situations in which this 
description applies is a separation of scales both in space and time. There are important 
cases however where local equilibrium apparently fails like aging phenomena in disordered 
systems due to insufficient ergodicity. 

The simplest non equilibrium states one can imagine are stationary states of systems 
in contact with different reservoirs and/or under the action of external fields. In such 
cases, to the contrary of equilibrium, there are currents (electrical, heat, matter of various 
chemical constitutions ...) through the system whose macroscopic behavior is encoded in 
transport coefficients like the diffusion coefficient, the conductivity or the mobility. The 
ideal would be to approach the study of these states starting from a microscopic model 
of atoms interacting with realistic forces and evolving with Newtonian dynamics. This is 
beyond the reach of present day mathematical tools and much simpler models have to be 
adopted in the reasonable hope that some essential features are adequately captured. 

In the last decades stochastic lattice gases have provided a very useful laboratory for 
studying properties of stationary non equilibrium states (SNS). Besides many interest- 
ing results specific to the different models considered, the following features of general 
significance have emerged. 

1. Local equilibrium and hydrodynamic equations have been derived rigorously from 
the microscopic dynamics for a wide class of stochastic models. 

2. A definition of non equilibrium thermodynamic functionals has emerged via a theory 
of dynamic large fluctuations, moreover a general equation which they have to satisfy 
has been established. This is a time independent Hamilton- Jacobi (H-J) equation 
whose independent arguments are the local thermodynamic variables and requires as 
input the transport coefficients. These coefficients can be either calculated explicitly 
for given models or obtained from measurements so that H-J can be used also as a 
phenomenological equation. 

3. Non equilibrium long range correlations, which have been observed experimentally in 
various types of fluids [22], appear to be generic consequences of H-J. An important 
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connection between the behavior of the mobihty and the sign of these correlations 
can be derived from the H-J equation. 

4. An analysis of the fluctuations of the currents averaged over long times has revealed 
the possibility of different dynamical regimes, which are interpreted as dynamical 
phase transitions. Such phase transitions have actually been proved to exist in some 
models. This theoretical prediction should be investigated experimentally. 

As an overall comment we may say that the macroscopic theory obtained so far encom- 
passes the theory developed long ago by Onsager [45] and then by Onsager-Machlup [46] 
for states close to equilibrium and we believe to be applicable in general to states where 
diffusion is the dominant dynamical mechanism. 

The present paper intends to provide a unified introduction to models studied inten- 
sively in the last decade for which several results have been obtained. We shall concentrate 
on time independent properties that is mainly on the topics [2] and [3] mentioned above. 
Our treatment is based on an approach developed by the authors in [2-5] . For item [T] we 
refer to [39,50,51] for the case of periodic boundary conditions and to [26,27] for open 
systems. For item|U which we do not discuss here, see [6-9,13-15]. For recent overviews 
on nonequilibrium phenomena see also [34,44]. 

There are two parts in the article according to the natural separation into microscopic 
and macroscopic properties. As a rule we do not include proofs except for statements 
which require a short argument. Most of the results discussed here are in published articles 
to which we shall refer for the details. We do outline however the following new results 
which will be the subject of forthcoming papers. In Section 12.41 we consider the KMP 
process [38]. We give an explicit representation of the invariant measure of this process in 
the case of a single oscillator and we compute exactly the two point correlations for the 
general case. In Section [3l3] we show, for a particular class of one-dimensional models, that 
nonequilibrium long range correlations are positive if the mobility is convex and negative 
if it is concave. The general case is discussed in [10]. In Section [3.41 we show that for 
any weakly asymmetric model with periodic boundary conditions the nonequilibrium free 
energy does not depend on the external field, so that it coincides with the equilibrium one 
and there are no long range correlations. In Section [3^ we consider the one-dimensional 
boundary driven totally asymmetric exclusion process. For a particular choice of the 
parameters, starting from the results in [23] we obtained a new variational representation 
of the nonequilibrium free energy. In particular, while the representation in [21] requires 
the maximization of a trial functional, we show how it can be formulated as a minimization 
problem. 

2 Nonreversible microscopic models 

Stochastic lattice gases are - loosely speaking - a collection of random walks moving 
in the lattice and interacting with each other. These "particles" are to be considered 
indistinguishable. Accordingly, the microscopic state is specified by giving the occupation 
number in each site of the lattice. The effect of the interaction is that the jump rates 
depend on the local configuration of the particles. For non-isolated systems we model the 
effect of the reservoirs by adding creation/annihilation of particles at the boundary. The 
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effect of an external field is modeled by perturbing the rates and giving a net drift toward 
a specified direction. 

In Section [2?T] we give the precise definition of non equilibrium stochastic lattice gases. 
Some special models, the zero range process and the exclusion process, are discussed in 
Sections 12.21 and 12.31 For the zero range process the invariant measure is always product 
and can be computed explicitly. On the other hand, the boundary driven exclusion process 
carries long range correlations, which can be computed explicitly in the one dimensional 
case. In Section 12.41 we recall the definition of the KMP process [38] and we compare it 
with the exclusion process. In Section [231 we consider gradient lattice gases with periodic 
boundary conditions; the peculiarity of such models is that the invariant measure does not 
depend on the applied external field. In Section 12.61 we consider the Glauber + Kawasaki 
model, in which a reaction term allowing creation/annihilation of particles in the bulk is 
added. We discuss under which conditions on the reaction rates it is reversible. Finally, in 
Section [2171 we consider the boundary driven totally asymmetric exclusion process and we 
recall the representation of the invariant measure obtained in [23]. This representation 
suggests a new variational expression for the nonequilibrium free energy that will be 
discussed in Section 13. 6[ 

2.1. Stochastic lattice gases 

As basic microscopic model we consider a stochastic lattice gas in a finite domain, 
with an external field, and either with periodic boundary conditions or with particle 
reservoirs at the boundary. The process can be informally described as follows. At each 
site, independently from the others, particles wait exponential times at the end of which 
one of them jumps to a neighboring site. In the case of particle reservoirs, in addition 
to this dynamics, we have creation and annihilation of particles, at exponential times, 
at the boundary. To define formally the microscopic dynamics, recall that a continuous 
time Markov chain Ut on some state space VL can be described in term of its infinitesimal 
generator L defined as follows. Let / : f2 ^ M be an observable, then 



where E( | ) is the conditional expectation, so that the expected infinitesimal increment of 
f{uJt) is {Lf){ujt) dt. The transition probability of the Markov process uot is then given by 
the kernel of the semigroup generated by L, i.e. 



Let A be the (i-dimensional torus of side length one, i.e. (M/Z)'^, respectively a smooth 
domain in M'^, and, given an integer A^ > 1, set A^r := (Z/A^Z)'', respectively A^r := 
(A^A) n U^. The configuration space is X^^, where X is a subset of N, e.g. X = {0, 1} 
when an exclusion principle is imposed and X = N when there is no limitation on the 
number of particles. The number of particles at the site a; G A^v is denoted hj rj^ & X 
and the whole configuration by ?7 G X^^. The microscopic dynamics is then specified 
by a continuous time Markov chain on the state space X^^ with infinitesimal generator 
given by Ltv = I/q.at, resp. Ljq = \Lq,n + -^^b.Af] , if A is the torus, resp. a smooth domain 



^f{u;t+H)\^t) = {Lf){uJt)h + o{h) 



(2.1) 
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in M'', where, for functions / : X^^ M, 

Lo,Nfiv) = I c^Av)[fi^''''v)-fiv)] (2.3) 

a;,yeAjv 

\x-y\=l 

Lb,Nm = ^ E {^'^M[fi^''''v)-fiv)]+CyAv)[fi^'^'v)-m]} (2.4) 

\x-y\=l 

Here |x| stands for the usual Euchdean norm. For x, y G An, a^^'^T] is the configuration 
obtained from 77 by moving a particle from x to y, i.e. 

772; - 1 if z = X 

and similarly, ifx G Ajv, y ^ An, then cr^'^?] is obtained from 77 by creating a particle at x, 
while (T^'^r/ is obtained by annihilating a particle at x. Therefore for x,y & An, Cx,y is the 
rate at which a particle at x jumps to y. We assume that Cx,y{rf) = if a^'^?] ^ X'^^ so 
that Lq^n and L^^at are well defined linear operators on the set of functions / : X'^^ R. 
The generator Lq^n describes the bulk dynamics which preserves the total number of 
particles whereas L^^at models the particle reservoirs at the boundary of Aat. 

We assume that the bulk rates Cx,y, x,y E An, are obtained starting from reversible 
rates c°, y satisfying the detailed balance with respect to a Gibbs measure defined by a 
Hamiltonian Ti, and perturbing them with an external field F . Likewise, in the case of 
particle reservoirs, we assume that the boundary rates c^^y, Cy^x, x E An, y ^ An, are 
obtained from rates c^y, c^^^ satisfying the local detailed balance with respect to H, and 
in presence of a chemical potential Aq, and again perturbed by the external field F. Our 
analysis is restricted to the high temperature phase, in particular we shall assume that 
the correlations in the Gibbs measure decay exponentially. 

The above conditions are met by the following formal definitions. Consider jump rates 
c° y satisfying the detailed balance with respect to the Gibbs measure associated to the 
Hamiltonian Ti : M with free boundary conditions. For the bulk rates this means 

cly{^) = exp { - [Hia^'yr^) - n{^)\ ] < ,(a^'''^) , x,y e An (2.5) 

Note that we included the inverse temperature in Ti. As before if o"^'^?7 ^ we assume 
'^xyiv) = 0. From a mathematical point of view, the detailed balance condition means 
that the generator is self-adjoint w.r.t. the Gibbs measure n{r]) oc e"^^^-*; namely if we let 
Lq ^ be the generator in (12.31) with c replaced by c°, for each /, g : X"^^ — >■ M we have 

(/, Ll^g), := J2 f'iv) fiv) L^NQiv) = (^o,iv/, 9), (2.6) 

V 

Let Atv := {x G Z"' I 3 y G An with \x — y\ < 1} be the 1-neighborhood of Aat. When 
Aat is the discrete torus we agree that Aat = Aa^. We also let OAn := An \ An- In the 
case when A is not the torus, the boundary dynamics with no external field is specified 
as follows. Denote by Aq : OAn K the chemical potential of the reservoirs. If x G Aat, 
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y ^ A]\r the detailed balance condition fl2.5p is modified by adding the chemical potential 
Ao: 

cl^iv) = exp { - [nia^'yr^) - Hiv)] - Ao(y) } < .(^^'^) , x e A^, y ^ Aj, (2.7) 

We denote by -B(Ajv) := {{x,y) \ x,y e Aj^ , {x,y} H A^ ^ ^ , \x — y\ = 1} the 
collections of ordered bonds intersecting A^r. A discrete vector field is is then defined as 
a real function F : B{An) —>■ M satisfying F{x, y) = —F{y, x) for any (x, y) G i3(Ajv). An 
asymmetric lattice gas is defined by the jump rates 

:= e^^^'y^ cl^iv) (2.8) 

where c° are the unperturbed rates and F is a discrete vector field. 
The case of weakly asymmetric models is obtained by choosing 

F{x,y) ^ F^{x,y) = e[^) ■ ^ (2.9) 

where E' : A — > M'^ is a smooth vector field and ■ denotes the inner product in M'^. Namely, 
for large, by expanding the exponential, particles at site x feel a drift N~^E{x/N). 

We can rewrite the full generator L^, using the notation introduced above, as follows 

(x,3/)GB(AAr) 

Fix an initial condition rj G X^^. The trajectory of the Markov process Tj{t)^ t > 0, is then 
an element on the path space D(M+;X^^), which consists of piecewise constant paths 
with values in X^^. We shall denote by the probability measure on 
corresponding to the distribution of the process ?7(t), t > 0, with initial condition rj. It is 
related to the generator L^v by (?7(t) = rj') = e^^^ {rj.rj'). 

A probability measure fi^ on X^'^ is an invariant measure for the process r]{t) if 

J2 Mv)e'''Hv,v')=Mv') (2.11) 

namely, if we distribute the initial condition rj according to /iat, then the distribution of 
r]{t) is fiN for any t > 0. According to general results on Markov processes, if the process 
is irreducible, i.e. there is a strictly positive probability to go from any state to any other, 
then the invariant measure is unique and it encodes the long time behavior of the system. 
More precisely, starting from any configuration rj the distribution of //(t) converges to fiN 
as t — > oo. In the case of a stochastic lattice gas with particle reservoirs, if Cx,y{rf) > 
for any rj G X"^^ and any (x, y) G B{An) such that a^'^rj G X"^^, then the process is 
irreducible and there exists a unique invariant measure. On the other hand, if A/^r is the 
discrete torus, the total number of particles YIix&Km ''^^ conserved and there exists a 
one parameter family of invariant measures. Since in general the transition probability 
cannot be expressed in a closed form, condition (12.111) is not convenient. However, it 
is easy to obtain a necessary and sufficient infinitesimal condition for a measure to be 



5 



invariant. The measure jj^N is invariant for the process generated by L^r if and only if, for 
any / : X^^ R, we have 

fiN{LNf)=0 (2.12) 

where hereafter for a measure /i and an observable / we denote by yu(/) the expectation 
of / with respect to /i. 

If the generator satisfies the detailed balance condition with respect to some mea- 
sure fiN, namely 

MgLNf) = MfLNg) Wf,g -.x^^ ^r (2.13) 

then /iTv is necessarily an invariant measure. In such a case the process is said to be 
reversible. This terminology is due to the following fact. Let P^^ the stationary process, 
i.e. the distribution on the path space induced by the Markov process with initial condition 
distributed according to the invariant measure fiN- Since /^at is invariant, the measure P^^ 
is invariant with respect to time shifts. We can thus regard P^ as a measure on paths 
defined also for t < 0, i.e. as a probability on D(R; X^'^). This probability is invariant 
under time reversal if and only if the measure /xat is reversible, i.e. 02.131) holds. More 
generally, if we denote by the time reversal, i.e. {'dr]) [t) := ri{—t), we have that P^^ o d 
is the stationary process with generator the adjoint to in L2{dfiN)- In particular, 
if ( ]2A3il holds, we have F^^ o ^ = P^^. 

When there is a unique invariant measure which is not reversible, we say the corre- 
sponding process is non-reversible. As we shall discuss, nonequilibrium models are nec- 
essarily non-reversible, while there exist non-reversible processes describing equilibrium 
phenomena, see e.g. [1,30-32]. The main topic that we shall discuss is the asymptotic 
behavior, as N diverges, of the invariant measure fiN for specific classes of non- reversible 
models. 

Conditions 02. 5p and 02.71) are called local detailed balance for the following reason. If 
the chemical potential Aq is constant, it is easy to show that the Gibbs measure fi{ri) oc 
^-'H(v)+XoJ2z'^^ ig reversible with respect to the generator L^, that is (12.101) with rates c^y. 
On the other hand, if Aq is not constant, the boundary dynamics forces a current in the 
system which becomes non reversible. 

We discuss next the effect of the external field. A very particular choice of F is that 
of a discrete gradient vector field, such that 

F{x,y) = ^-[X{y)-X{x)] (2.14) 

for some function A : A^r — » M. If we further assume that \{y) = \o{y) for y G dA^, recall 
that Aq is the chemical potential of the boundary reservoirs, then it is easily shown that 
the generator L^r in O2.10p is reversible w.r.t. the measure 

i^Niv) = 4r exp I - niv) + Yl ^(^) ^4 (2.15) 



where is the appropriate normalization constant. In this situation the reversibility of 
the process is due to the fact that the driving from the reservoirs and the external field 
compensate. 
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An equivalence principle. We can interpret the above result from two different perspectives, 
getting the answer to two opposite questions. 

Consider a model with a given chemical potential Aq at the boundary. We ask if we can 
find an external field F which compensates the driving from the boundary, namely such 
that the corresponding stationary measure /xat is reversible. The answer to this question 
is certainly yes. In fact, from the above consideration, we have a whole family of external 
fields that fulfill this condition: take a gradient vector field F as in fl2.14l) with A such that 
X{y) = Xo{y) for y G dAj^. In this case the corresponding stationary measure happens to 
be the Gibbs measure (I2.15p . 

Conversely, suppose that we have an asymmetric model with a given external field 
F. We ask if we can find a chemical potential Aq : dA^ — > M such that the model is 
reversible. We can immediately answer affirmatively this question if we know that the 
external field F is gradient, i.e. (12.141) holds for some A : A^r — > M. In this case we can 
just fix, up to an overall additive constant, Xo{y) = X{y) , y E dApf. 

2.2. Zero range process 

The so-called zero range process is a special case of the lattice gases introduced in 
Section 12. 1[ In each site any number of particles is allowed so that X = N and the bulk 
symmetric jump rates are 

cl,y{v) = giVx) , x,y eAN,\x-y\ = l (2.16) 

where g : N ^ M_|_ is a function satisfying g{0) = with at most linear growth. In other 
words, the jump rate from x to y depends only on the occupation number at x; this 
explains the name of the model. Given Aq : dAj^j — > M, we choose the boundary rates as 

clyiv) = 9iv.) , c°,(r/) = e"«(^) , x G A^, y G M^, |a; - y| = 1 (2.17) 

It is not hard to check that the detailed balance conditions (12.51) and (12.71) are satisfied 
with the Hamiltonian 

where, by definition, g{0)\ = 1 and g{k)\ = g{l) ■ ■ ■ g{k) for k>\. The particular case in 
which g is the identity, i.e. g{k) = fc, G N, corresponds to independent random walks 
described in terms of the occupation variables rj G N^^. 

A peculiar feature of this model is that its invariant measure, both with particle 
reservoirs and external field, is always product. Consider this model with external field 
F. Denote by : A^v IR+ the solution to 



E 



[e^(^'"V(y) - e^("'^V(a;)] = , a; G A 



N 



y&^N (2.18) 

\x-v\=l ^ ' 

•0(0;) = exp{Ao(x)} , x E dA^ 

The invariant measure of the zero range process with external field F and boundary 
chemical potential Aq is then the grand-canonical product measure = HxeAjv /^^ with 
marginal distributions 

1 ihix)^ 

^^AVx = k) = — — -- ^ (2.19) 
Z(^(x)) ^(A;)! 
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where 

oo ^ 

^(^) = g^ (2^20) 

is the normahzing constant. This can be verified by showing that (12.121) holds. If A^v is the 
discrete torus and F vanishes, any constant ip solves (12.181) . the corresponding invariant 
measures are thus the grand-canonical measures with arbitrary chemical potential. Since 
the invariant measure is always product, the zero range process never exhibits long range 
correlations. 

2. 3. Exclusion process 

The exclusion process is a much studied stochastic lattice gas. In this model an 
exclusion principle is imposed. In each site x G A^r at most one particle is allowed so that 
X = {0, 1} and there is no other interaction. The symmetrical bulk rates are defined by 

c°,3/(^) = Vx{'^ -Vy), e An, \x - y\ = 1 (2.21) 

namely a particle at x jumps to a nearest neighbor site y with rate 1/2 if that site is 
empty. Then (12.51) holds with H = 0. Note that the rates (12.211) satisfy the constraint 
c^^y{r]) = if a^'y-T] ^ {0,1}^'^. Given a chemical potential Aq : dA^ R, the local 
detailed balance condition (12.71) is met by choosing the boundary rates as 

clyiv) = vJ<iy) , cl^iv) = (1 - V.)Kiy)e^'>^y^ , x e A^, |/ e OAn, \x-y\ = l (2.22) 

for some K : dA^ — > M+. 

We first discuss this symmetric case. In the case of periodic boundary conditions 
there is a one parameter family of invariant measures which are the Bernoulli measures 
with an arbitrary parameter. Since the total number of particles is conserved, given 
A; > 0, we can consider the process on the set T.^^^ = {v ^ {0; l}'^^ : 12x&An ~ ^J"' 
In this set the process is irreducible and the unique invariant measure is the uniform 
measure on T^^^k which is the canonical ensemble associated to the Bernoulli measures. 
In the case with particles reservoirs, if the chemical potential Aq is constant then the 
unique invariant measure is the Bernoulli measure with parameter p = e'*'"/(l + e'^°), i.e. 
jJ'Niji) = e^^^^ehf^'nx _|_ gAo-jiAivi^ both these situations the process is reversible. 

One dimensional boundary driven exclusion process. Unlike the zero range model, if Aq 
is not constant, so that this becomes a nonequilibrium model, the invariant measure is 
not a product measure and carries long range correlations. Let us discuss in more detail 
the one-dimensional case. Assume that A = (0, 1) so that A^v = {1,...,A^ — 1}; we also 
let Aq := Ao(0) and Ai := Ao(iV) be the two chemical potentials of the reservoirs. 

An old result by Kingman [37] computes the marginals of the unique invariant measure 
for a special choice of the injection rates. More precisely, in the case analyzed by Kingman 
the bulk rates are as in (I2.2ip while the boundary rates are obtained by the following 
hmiting procedure. In flZ^ choose K{Q) = (e^ + e-^)-^e-^o/2 i^(Ar) = (e^ + 
e~"^)~^e~^^^'^ for some A G M. Consider then the asymmetric model with rates c.j.y 
as in (12. 8p by introducing the external field F given by -F(0, 1) = F{N — 1, A^) = A, 
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F(1,0) = F{N,N — 1) = —A, and F{x,y) = in all the remaining bonds. Finally we 
take the limit A ^ oo obtaining 

CoAv) = (1 - Vi) e^"/' , CM^i.Niv) = VN^i e-^^/' , = cn,n-i{v) = (2.23) 

i.e. from the left endpoint particles enter with rate (l/2)e^°^'^ but do not exit, while 
particles from the right endpoint exit with rate (l/2)e~'^^/^ but do not enter. 

By some smart duality computations, Kingman shows that, for this particular choice 
of the boundary rates, the marginals of the invariant measure fiN are 

, . {A - m - xi){A - m - X2 + I) ■ ■ ■ {A - 1 - x.m) „ 

iVx^ = 1, ■ ■ ■ , r^x„ = 1 = 7^ TT^ — T 7^ — -T 2.24 

[B — m)[B — m + 1) ■ ■ ■ [B — I) 

where 1 < a;i < X2 < ■ ■ ■ < Xm < N — 1 are lattice sites and the parameters A and B are 
defined as 

A = N + e^i/^ ; B = N -1 + e^^/^ + g-Ao/2 ^2.25) 

More recent work based on matrix methods, allows to get some representation of the 
invariant measure in the general one- dimensional case, see e.g. [42, 48] and references 
therein. 

We consider now the one-dimensional boundary driven symmetric exclusion model 
with boundary rates as in ^TM with K{<d) = (1 + 6^°)"^ and K{N) = (1 + e^^)'^. 
As before Aq and Ai are the chemical potentials of the boundary reservoirs. Letting 
Pi = e^^ /{I + e^^), i = 0,1, be the corresponding densities, we then get 

Clfiiv) = (l-Po)??! , Co,l(?7) = po(l-'7l) , C7V_i,Ar(?7) = (l-pi)?77V_i , Cn,N-1 = Pl(l-?7Ar-l) 

Let /iAT be the unique invariant measure, it is not difficult to show that the density 
profile Hn{Vx) is hnear so that 

X 

fiNiVx) = Po + ^ (Pi - Po) (2.26) 

As first shown in [49], it is also possible to obtain a closed expression for the two-point 
correlations. For l<x<?/<A^ — Iwe have 

To prove this result it is enough to compute LN^rjxrjy), i.e. the action of Lat on the function 
rixTiy, and solve the equation 

PN{LNiVxVy)) = (2.28) 

Note that, if we take x < y at distance 0{N) from the boundary, then the covariance 
between rj^ and riy is of order 0{1/N). Moreover the random variables rjx and rjy are 
negatively correlated. This is the same qualitative behavior of the two-point correlation 
for the uniform measure on Sat^. As we shall show below, quite the opposite behavior is 
found in another model, the KMP process. 
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One dimensional periodic asymmetric exclusion process. We finally discuss the case of 
the asymmetric exclusion process on the discrete torus. It is defined by the jump rates 



Cx,x+i{v) = e Vxi^ - Vx+i) c^+iAv) = e Vx+i{^ - Vx) (2.29) 

for some F G M so that (12.81) holds with constant external field F. A simple computation 
shows that the Bernoulli measure fip with arbitrary p G [0, 1] is an invariant measure. 
Note however that for F the process is not reversible; in fact the stationary process 
w.r.t. /ip carries the mean current p(l — p) sinh(F). Unlike the zero range, if the external 
field F is not constant the invariant measures are in general not anymore product. Note 
however that if F is a gradient vector field, as shown before, the process is reversible 
w.r.t. a product measure. We emphasize that a constant vector field on the torus is not 
gradient. 

2.4- The boundary driven KMP process 

The Kipnis-Marchioro-Presutti (KMP) model [38] describes a chain of one-dimensional 
harmonic oscillators which are mechanically uncoupled but interact stochastically as fol- 
lows. Each pair of nearest neighbors oscillators waits an exponential time of rate one and 
then redistributes uniformly its total energy. The two oscillators at the end points are 
coupled to heat reservoirs. Since the single spin space state is not discrete and the elemen- 
tary dynamics is associated to the bonds, this model does not really fit in the framework 
introduced in Section [21 However the precise definition of the model is straightforward. 
Let A = (0, 1) so that Aat = iVA n Z = {1, . . . , - 1}. We denote by the energy of 
the oscillator at the site x G Aat, so that the state space is M.^^ . On it we introduce the 
Markov generator L^r as follows. Given (x, y) G B{Kn) and p G [0, 1] we let ^(^'^^^'^ be the 
configuration obtained from ^ by moving a fraction p of the total energy + across 
the bond (x, y) to x and a fraction 1 — p to y, i.e. 

{iz ii x,y 

p{^x + Q if z = x 

{l-p){^. + Q if z = y 

We then set Ljy := X^^o^ where, for / : R^^^ — R, the bulk dynamics is given by 

Lx,.+J{0 ■■= fdp [/(e(-'-+i)'^) - /(O] , X = 1, . . . , AT - 2 
while the boundary generators Lq i and L^_i^^ are 

Lo,if{0 ■■= / rfeo;^e"^ / rfp[/(e(°'^)'^)-/(0] 
Jo ^0 Jo 

L^-i,Nf{i) := / rfe^;^e-^ / rfp[/(e(^-^'^)'^)-/(0] (2.30) 
Jo J-i Jo 

Namely, we suppose that there is an energy exchange across the ghost bonds (0, 1) and 
(A^ — 1,A^), and we put at the sites and A^ oscillators whose energies are randomly 
chosen according to the Gibbs distributions with temperatures Tq and Ti. 
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We emphasize that the above choice of the boundary dynamics differs shghtly from the 
original one in [38]. Besides being more natural, this choice simplifies some microscopic 
computations. Note that in the case T = Tq = Ti, namely of an equilibrium model, the 
above process is reversible with respect to the Gibbs measure 

c?/^^(0= n T^"^^^"^^^- (2.31) 

x=l 

which is just the product of exponential distributions. 

Later, in order to find a closed expression for the microscopic two-point correlation 
functions, we introduce a more general class of boundary dynamics which is obtained by 
replacing in f l2.30p the two exponential distribution on boundary sites and by other 
two probability measures on M4.. Of course the macroscopic behavior is the same for any 
reasonable choice of the boundary dynamics. 

Invariant measure for a single oscillator. We consider the KMP model with a single 
oscillator, i.e. N = 2. Even in this the system is in thermal contact with two 

reservoirs, its stationary state is not trivial. We next show that the invariant measure is 
a mixture of the Gibbs distributions with temperatures between Tq and Ti. Furthermore 
we compute the weight of each distribution which turns out to be the arcsin law in the 
interval [To,Ti], here we assume Tq < Ti. We emphasize that this result depends on the 
specific choice of the boundary dynamics. 

We claim that the invariant measure (a probability measure on M4.) is absolutely 
continuous w.r.t. the Lebesgue measure and its density can be expressed as 

dfi f^' , 1 



I ^dQT,,TAT)^e~^'^ (2.32) 



where ^To,Ti is the arcsin distribution in the interval [Tq, Ti], namely for T in this interval 
we have 

dQT, Ti (T) = - dT (2.33) 

' ^ v/(ri-T)(T-To) 

To show that /i in f l2.32p is the invariant measure of the KMP process with a single 
oscillator, we need to check that for each smooth real function / on M+ (12.121) holds. By 
linearity and approximation by linear combinations of exponential functions, it is enough 
to show that (12.121) holds if /(^) = exp{— A^}, A > 0. With this choice we have 



^0 JQ Ti 







dp 







-e 



Ll + ApTo 1 + ApTi 



If we now take the average of the above expression when ^ is an exponential random 
variable of parameter T we get 



'fe-^/^L/(0= / dp 



^ ^ 1 1 







:i + ApTo)(l + ApT) {I + XpT){l + XpT^) 1 + ATJ 

(2.34) 
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We next note that the arcsin distribution in the interval [Tq, Ti] is characterized by 
the following property. For each 7 > we have 



dgTo,Ti{T) 



I + 7T ^(l + 7To)(l + 7Ti 



(2.35) 



The integral on the l.h.s can be in fact computed by using the density in fl2.33p and the 
residue theorem. Conversely, by expanding the above equation in power series of 7, we 
get that the moments of ^Tq.Ti are determined. 

Recalling fl2.34p . to complete the proof of (12.321) it remains to show that 



d^ToMT) / dp 
Jo 



1 



1 



;i + ApTo)(l + XpT) (1 + ApT)(l + XpT^] 
which in view of (I2.35P is equivalent to 

1 1 



1 + AT 







dp 



+ 



:i + Ai^ro)3/2(l + ApTi)V2 (l + ApTo)i/2(i + Apri)3/2_ 
2 



(l + ATo)V2(i + ATi)V2 
By a direct integration we get 

dp ■ 



1 + ApTo)3/2(l + XpT,y/2 A(Ti - To) 



yi + ATi 

yi + ATo 



1 



and simple algebraic computations yield the result. 

It seems quite hard to obtain an analogous representation for > 3. On the other 
hand we compute explicitly below the one and two point correlation functions of /^tv for 
any N > 2. 

Two point correlations. We here consider the KMP process with boundary dynamics given 
by 

POO p1 

Lo,if{0 ■■= / rf^o^(eo) / cip [/(e(°'^)'^) - /(O] 
Jo Jo 



L 



N-1,N 



/(O := 



di^^i^N) / dp [/(e 



■:iN~l,N),p^ 



where uj^ , i = 0,1, are probability measures on IR+ with mean Tj and variance 



,N 



(Ti - To 

N{N + -. 



0,1 



(2.36) 



and note that the exponential distributions chosen in (I2.30p fail to satisfy the above 
condition only by a term 0{N~'^). 

Given Tq < Ti, let he the invariant measure of the KMP process with — 1 
oscillators and set E^{x) := fiNi^x), x = 1, ■ ■ ■ , N — 1, as well as En{0) := Tq, Eiy{N) : = 
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Ti. By choosing linear functions / in fl2.12p and computing Lj^^^, x = 1, ■ ■ ■ , N — 1, we 
get a closed equation for E]\r which yields 

E^(x)=To + (Ti-To)^ (2.37) 

Let Cn{x, y) := fiNi^x] Q = /i7v(^x ^y) - EN{x)EN{y), x, y G {1, ■ ■ • , iV - 1}, be the 
two point correlation function of yU AT. We also set CAr(0, 0) := {\^ — Tq\^^, Cn{N,N) : = 
- ^i]^), CN{0,y) = Cn{x,N) := foi 1 < y < N, < X < N - I. By choosing 
quadratic functions in fl2.12p . by some elementary but tedious computations we get that 
Cn{-, ■) solves 

' {A^ + A^)CNix,y) = l<x<y<N-l, y-x>2 

Cn{x, X + 2) + Cn{x - 1, X + 1) - ^ Cn{x, X + 1) + I Cn{x, x) 

+ |C;v(x + l,x + l) = I (^i^)' + |E^(x)Ejv(x + l) l<x<A^-2 

Cn{x - 1, X - 1) + CAr(x + 1, X + 1) + 2 Cn{x - 1, x) + 2 Cn{x, X + 1) 

-4Cn{x,x) = -2{^^Y -2En{x)^ 1<x<N-1 

where f{x,y) = /(x + l,y) + /(x — l,y) — 2/(x, y) is the discrete Laplacian w.r.t. x 
and is the discrete Laplacian w.r.t. y. 

As can be easily checked, the solution is given by 

< X < y < N 

(2.38) 

) + ^WtM 0<x = y<iV 

Comparing fl2.38p with (12.271) we observe that the off diagonal terms are essentially 
the same in the macroscopic limit N ^ oo. We emphasize however that the sign is 
different: while for the boundary driven symmetric exclusion the occupation variables rj^, 
X e Aat, are negatively correlated, for the KMP process the local energies C,^, x G A^v, 
are positively correlated. As we shall discuss in Section 13.31 this qualitative difference is 
related to the different convexity properties of the mobilities of the two models. For the 
exclusion it is concave while it is convex for KMP. 

We mention that an analogous computation has been recently performed for a some- 
what similar model, see [33]. 

2. 5. Gradient models with periodic boundary conditions 

In this section we consider only the case of periodic boundary conditions, namely A^v 
is the discrete torus (Z/A^Z)'^. We also assume that the model is translationally covariant 
in the sense that, for any (x, y) G B{Aiy), z G Atv, and r] G X^^, we have 

= c^+z,y+z{rzV) (2-39) 
where is the space shift, i.e. {Tzrj)x '■= rj^-z- 



CN{x,y) 



{Ti-Tof X 
N+l N 



(1 



N 



En(x)^ + 2 



N+l N 



(1 
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Let the bulk rates c^y satisfy fl2.5l) . The expected instantaneous current across the 
bond {x,y) G B{An) is, up to a factor 2, 

The corresponding lattice gas (with no external field) satisfies the gradient condition if 
the discrete vector field j^yij]) is gradient for any 77 G X^^, namely there exist functions 
: ^ M, X G Aat, such that for any {x,y) G B{Kn) 

Jlyiv) = hyiv) - hM (2.40) 

The zero range process of Section [2l2] is a gradient lattice gas for any choice of the function 
g. Indeed, for the rates fl2.16l) condition fl2.40p holds with hx{ri) = —g{ri^). Also the 
exclusion process is gradient, fl2.40p holding with hx{ri) = —rj^. 

In the stochastic gases literature, see [39,50], the gradient condition is usually stated 
in a stronger form. More precisely, one considers a translationally invariant lattice gas 
on the whole lattice Z*^ and says that the model is gradient if there exists a function 
h : M which is local, i.e. it depends on r/^. only for a finite number of x G Z'^, and 

such that for any [x, y) G BilJ^) and 77 G X^'* 

3ly{ri) = KryV) - Kr.v) (2-41) 

Of course fl2.4ip implies fl2.40p for N large enough. Conversely, it is possible to show that 
if f l2.40p holds then there exists a function h : X^'^ —>■ M such that (I2.4ip holds for any 
(x,y) gS(A^). 

Consider now a lattice gas with constant (non zero) external field F. By this we mean 
that F{x,x±ei) = ±Fj, where Cj, i = 1, ■ ■ ■ , c? is the canonical basis in R*^ and {Fi, . . . , Fd) 
is a vector in M.'^. As discussed in [36], if the bulk rates c° satisfy the gradient condition 
fl2.4ip then the grand-canonical Gibbs measures exp{— 7^(77) + A ^^.g^^^ rj^}, A G M, which 
are the invariant measures for the system with no external field, are invariant also for 
the process with external field F, i.e. with rates Cx,x±e, = ^_^^.e^^\ i = 1, ■ ■ ■ , d. In 
particular this result shows that gradient lattice gases with constant external field and 
periodic boundary conditions do not exhibit long range correlations. In Section 13.11 we 
show that, from a macroscopic point of view, any weakly asymmetric lattice gas with 
periodic boundary conditions does not have long range correlations. 

We next discuss, from a microscopic point of view, gradient lattice gases in some 
more detail obtaining the above mentioned result as a particular case. Let us consider 
an asymmetric lattice gas with external field F : B{An) — > R, rates c^^y as in (12.81) . and 
generator given by (12.101) . We look for an invariant measure of the form (I2.15P for some 
A : A^r M. The condition for a stationary state is 

^'Niv) E c,,y{v)[f{^''''v)-f{v)]=0, \/f:X^-^R. (2.42) 

Performing some change of variables and using the conditions (12.50 and (12.71) of local 
detailed balance, (I2.42p becomes 

E f(v)f^Niv) E e-^(^)e^(^'^) [e'^'^clM - e'^^^cl^iv)] = (2.43) 

r]&X^N (a:,j/)GB{Ajv) 
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Let 

Note that if F is a discrete vector field, i.e. it satisfies F{x,y) = —F{y,x), then is not 
a discrete vector field but G'*' and are indeed discrete vector fields. We then get that 
( I2.43P is equivalent to 

E G\x,y)j^^M = 0. Vr/eX^-. (2.44) 

{x,y)eB{Ajv) 

Notice that f l2.44p is an orthogonality condition. In general there is no solution to fl2.44p : 
note in fact that it is a system of IXI'"^^' equations (corresponding to different particles 
configurations) but we have only |Ajv| parameters (corresponding to the chemical potential 
profile A : Ajv ^ K). Non-existence of solutions to (12.441) means that the invariant measure 
is not of the form f l2.15p . There are however few remarkable cases in which (12.441) can be 
easily solved. 

If the model is gradient, so that (I2.40p holds, we claim that fi^irj) in (I2.15P with A G M 
constant is an invariant measure for any asymmetric lattice gas provided the external field 
F satisfies 

Je^(--3/) _ e^(s''-)] =0 WxeAn (2.45) 

y:\x-y\=l 

that is the discrete vector field G^, A G M, has vanishing discrete divergence. Conversely, 
if we require that fi^{ri), A G M, is an invariant measure for any external field satisfying 
(I2.45p . we get that the rates c° ^ have to satisfy (I2.40p for some functions '■ X^^ M, 
X G Ajv. The proof of both statements are accomplished by some computations which 
essentially amounts to prove the Hodge theorem in a discrete setting, see [43]. 

Generalized gradient models. Consider asymmetric lattice gases with constant external 
fields F. Some computations show that fi^ as in (I2.15P is an invariant measure for any 
constant A G M if and only if the rates c° y satisfy 

E^'°-+e.(^)=0, Vr/GX^-,Vz = l,...,d. (2.46) 

which is exactly the condition that identifies the orthogonal complement, w.r.t. the inner 
product defined in (I2.44p . of the constant vector fields. Moreover (12.460 is equivalent 
to the following generalized gradient condition. There are function hij : X"^^ M, 
i, j = 1, ■ ■ ■ ,d, such that 

d 

jlx+eS^) = Yl [^i,ji^x+ejV) - hi,j{T^r])] , z = l,--- ,d (2.47) 

We finally mention that (12.471) is a particular case of the condition stated in [39, Def. 2.5]. 

To summarize the previous discussion, gradient models in the sense of (I2.40p have 
the property that any external field satisfying (12.450 will not change the invariant mea- 
sure, while generalized gradient models in the sense of (I2.47P have this property only for 
constant external fields. 
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2.6. Glauber + Kawasaki dynamics 

Unlike the models discussed so far, the so-called Glauber + Kawasaki process is not a 
lattice gas in the sense that the number of particles is not locally conserved. A reaction 
term allowing creation/annihilation of particles is added in the bulk. We consider the 
case with exclusion rule so that X = {0,1} and discuss only the one-dimensional case 
with periodic boundary condition, A^v a ring with sites. The generator is defined as 

LnM = \ E E ^M[f\^''v)-f{v)\ (2.48) 

(a;,3/)eB(Ajv) a;eAiv 

where denotes the particle flip at x, i.e. (cr^r/)^. = 1 — r]^ and {(J^ri)y = rjy for y ^ x. 
The first term of the generator corresponds to the symmetric exclusion process while the 
second one involves the reaction defined by the corresponding rates Cx, x G Ajy. The factor 
A^^ in (12.481) has been inserted to get, after diffusive rescaling, a meaningful macroscopic 
evolution. 

The first question one can ask is when there exists a reversible measure for this process. 
As we shall see, this happens only if we impose some restrictions on the reaction rates 
Cx- The condition of reversibility w.r.t. the measure is (12.61) . which in this case, after 
some algebra, reads 



E E air]) /(^"^) [cMMv) - Cx{a^v)M^'v)] = (2.49) 

where r/^'^+^ denotes the configuration obtained from t] by exchanging the occupation 
numbers in x and x + 1. Since this equality must hold for every g and /, this condition 
is equivalent to 

r /i7v(^) - /i7v(^'^'^+^) = 

\ Cx{r])lJ^N{v) - CxifT'^v) I^Nicr'^'n) = 
for any rj and x. The first condition imposes that the measure has the form 

/i^(r/) = M^( J2 V.) (2.51) 

namely iin must assign an equal weight to configurations with the same number of par- 
ticles. The second condition, with a /i at of this type, is a restriction on the reaction rates 
and on the function M^- The most general form of Cx{r]) that satisfies this condition is 

Cxiv) = Ai{l- rjx) h{Txr]) + A VxhirxV) (2.52) 

where Ai,A2 are arbitrary positive constants, and h : {0,1}'^^ is an arbitrary 

positive function such that h^a^rf) = hirf), i.e. it does not depend on 779. Recall that Tx 
denotes the shift by x. Notice that the rates Cx{ri) in (I2.52p are translation invariant, 
namely they satisfy Cx{rj) = cq^TxT]). With this choice, the unique reversible measure fiN 
is the Bernoulli measure with parameter p = ^^+3^' ["^2]. 
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We emphasize that periodic boundary conditions are crucial for the vahdity of (12.521) 
with a nontrivial h. In this special case there are no long range correlations. In Section 
we show that if the rates Cx are not of type (I2.52p then - generically - there are long range 
correlations. 

2. 7. Totally asymmetric exclusion process 

The one dimensional totally asymmetric exclusion process is the particular case of the 
one-dimensional asymmetric exclusion process introduced in Section in which particles 
jump only to the right. As discussed there, in the case of periodic boundary conditions, 
the invariant measures are the Bernoulli measures with any density. We instead consider 
here the boundary driven model. As usual we set Aa? = {1, ■ ■ ■ , — 1} and we let Aq and 
Ai be the chemical potentials of the two reservoirs. The bulk jump rates are 

Cx,x+i{'n) =Vx{'^-rix+i) , c^+iA'n) = ^ x = l,---,A^-2 (2.53) 

while the boundary rates are 

Cq,i{v) = r/ie 2 , cn-i,n{v) = VNe~~ , Ci,o(?7) = cn,n~i{v) = (2.54) 

These rates can be obtained from our standard choice by a limiting procedure analogous 
to the one described to get (I2.23p . 

The unique invariant measure for this model has an interesting representation due to 
Duchi and Schaeffer [23] that we briefly recall. We duplicate the variables by introducing 
new random variables ^ G {0, l}^'^. We then define a joint distribution for the variables 
(?7, ^ ) as follows. Let 

X 

E,:=J2iVz + Q-x, x = l,---,N-l 

z = l 

and Eq := 0. The measure uj^j gives positive weight only to complete configurations, 
defined by the conditions 

En-i=0, E^>0, x = 1,--- ,N -1 (2.55) 

Given a complete configuration we give some labels to the lattice sites according to 
the following rules: 

X G Atv has label W if C,x = and -Ex-i = E^ = 0; 

X G Atv has label 5 if = 1, E^-i = and there are no sites on the left of x labeled 
W. 

Let us denote by N^y = N]Y{rj,^) the number of sites with label W for the complete 
configuration (rj,^) and by Nb = NB{rj,^) the number of sites with label B. The measure 
pn is then defined as 

MV, = ^ exp \Nw Ai/2 - Nb X0/2} (2.56) 
where = Zi\f{Xi, A2) is the appropriate normalization constant. 
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The invariant measure of the boundary driven totally asymmetric exclusion process is 
then the first marginal of the measure u^, i.e. 

Mv)= Yl ^^(^'^) (2.57) 
5e{o,i}Aiv 

This result is proven by constructing a suitable Markov dynamics on the complete config- 
urations (?7, such that its projection to the t] variables coincides with the dynamics of 
the totally asymmetric exclusion process. The invariant measure of the enlarged Markov 
dynamics can be easily computed and yields f l2.56p . 

3 Macroscopic theory 

As previously stated, an issue that we want to discuss is the asymptotic behavior of the 
invariant measure fijsf. Let us first briefly recall the situation of reversible models. For 
definiteness consider a stochastic lattice gas with reservoirs at the boundary and assume 
that the chemical potential Aq of the boundary reservoirs is constant and that there is 
no external field. As discussed in Section 12.11 the unique invariant measure is the grand- 
canonical Gibbs distribution 

^^N{v) = - n{v) + Xo Yl ^4 (3-1) 

where, letting S^v,*: := {v G X^^ | J^xgAn ~ grand-canonical partition function 

Zn{X) is 

We then define po{X) as 

Po(A) := lim -^log/i^'fe^^-eA^ ''") (3.3) 

where lA^vl is the number of sites in A^v- Note that po can be easily related to the 
pressure. Let in fact po{X) := lim7v->oo IAat]"^ log Zjv(A) be the pressure, then po{X) = 

PoiXo + A) -po(Ao)- 

We then define the free energy /o as the Legendre transform of po, i.e. 

/o(p) :=sup{Ap-po(A)} (3.4) 

Asm 

According to the normalization chosen /o is a convex function which takes its minimum at 
the density associated to the chemical potential Aq, i.e. at po = Po(0) — P'oiXo)- Moreover 
/o(po) = 0. 

According to the Einstein fluctuation formula [24,40], see also Lanford's lectures [41] 
for a complete mathematical treatment, the free energy /q gives the asymptotic probability 
of observing a fluctuation of the density, namely 

^'nIt^ E ^-^^) ~exp{-|Ajv|/o(p)} (3.5) 
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here a ^ b means closeness in M and ~ denotes logarithmic equivalence as lA^vl diverges. 

In discussing nonequilibrium models, which are not translationally invariant, it is 
important to establish a generalization of the above fluctuation formula. We want to 
compute the asymptotic probability of a fluctuation not of the average density but of 
the density profile. In fact, already Einstein [24], considered density profiles in small 
fluctuations from equilibrium. We introduce the empirical density as follows. To each 
microscopic configuration r] G X^^ we associate a macroscopic profile n'^{u) = vr^(?7; u), 
n G A, by requiring that for each smooth function G : A ^ R 

(tt^, G)= jdu 7r^(n) G{u) = ^ G{x/N)r^, (3.6) 

so that 7T^{u) is the local density at the macroscopic point u = x/N in A. Let p = p{u) 
be a given density profile. Then fl3.5p can be recast as 

/^^"K ^ p) ~ exp { - N'^Mp)} (3.7) 

Here p ^ p' means that their averages over macroscopically small neighborhoods are close 
and J-'o{p) is the local and convex functional 

^o(p)= [dufo{p{u)) (3.8) 

J A 



For non-reversible systems we shall look for a fluctuation formula like (13.71) which, in 
the same spirit as Einstein, we shall consider as the definition of the nonequilibrium free 
energy. While in the reversible setting discussed above the invariant measure //at is given 
by the Gibbs distribution (13. ip . in a non-reversible system p^ is not, in general, explicitly 
known. For special models, powerful combinatorial methods have been used [19-21,25]. 
In the sequel we shall discuss instead the strategy introduced in [2,3] which is based on the 
following idea. As diverges the evolution of the thermodynamic variables is described 
by a closed macroscopic evolution called hydrodynamic equation. The microscopic details 
are then encoded in the transport coefficients appearing in the hydrodynamic equation. 
In the cases discussed here, these transport coefficients are the diffusion coefficient and the 
mobility. For the Glauber + Kawasaki dynamics also the reaction rates are involved. We 
then compute the asymptotic probability of fluctuations from the typical hydrodynamical 
behavior generalizing to a dynamical setting the Einstein fluctuation formula (13.70 . The 
nonequilibrium free energy JF is then characterized as the solution of a variational problem, 
from which we derive a Hamilton- Jacobi equation involving the transport coefficients. 
This is an infinite-dimensional strategy analogous to the Freidlin-Wentzell theory for 
diffusion processes, [29]. 

Of course, in the case of reversible systems, the solution to the Hamilton- Jacobi equa- 
tion coincides with the equilibrium free energy J-'q. This is essentially the characterization 
of JFq given by Onsager-Machlup [46], extended to a non-linear context. 

In Section [3TT] we discuss the hydrodynamics and the associated dynamical large devi- 
ations principle of weakly asymmetric lattice gases. In Section [32] we recall the derivation 
of the Hamilton- Jacobi equation and we discuss the form of the nonequilibrium free en- 
ergy for the specific models introduced in Section [21 We will also discuss a toy model for 
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the invariant measure of the KMP process. In Section 13.31 we obtain macroscopic equa- 
tion satisfied by the correlation functions and we discuss whether correlations are positive 
or negative. In Section 13.41 we show that for weakly asymmetric lattice gases with peri- 
odic boundary conditions the nonequilibrium free energy coincides with the equilibrium 
one. In Section 13.51 we discuss the macroscopic property of the Glauber + Kawasaki 
dynamics [1,32]. Finally, in Section [XUl starting from the results in [23], we show how 
the representation for the nonequilibrium free energy of the totally asymmetric exclusion 
process obtained in [21] can be formulated as a minimization problem. 

3.1. Hydrodynamics and dynamical large deviations 

We consider an asymmetric model as defined by the rates (12.81) . If the microscopic 
external field F is of order 1, the appropriate scaling is the Euler one, i.e. both space 
and time are rescaled by a factor A^, and the hydrodynamic equation is given by an 
hyperbolic equation, see [39] and references therein. We here consider instead the case 
in which the external field is of the order 1/N as in fl2.9p . Then the hydrodynamic limit 
is obtained in the diffusive scaling and given by a parabolic equation. Let vr^(t) be the 
empirical density, as defined in (13. 6p . corresponding to the particles configuration at time 
A^^t, 7r^(t, m) is then a random space-time trajectory; as — > oo it converges however 
to a deterministic function. Referring to [39,50,51] for periodic boundary conditions and 
to [3,7,26,27] for open systems, we here state the law of large numbers, as — oo, of 
the empirical density tt^ for weakly asymmetric lattice gases. The macroscopic evolution 
of the density is described by a (in general nonlinear) diffusion equation with a transport 
term corresponding to the external field, namely 



where D is the diffusion matrix, obtained from the microscopic dynamics by a Green- 
Kubo formula [50, II. 2. 2], and x is the mobility matrix, obtained by linear response 
theory [50, II. 2. 5]. In (13. 9p ■ denotes the standard inner product in M'^. This equation has 
to be supplemented by the boundary conditions which are either periodic when A is the 
torus or the non-homogeneous Dirichlet condition 



in the case of boundary driven systems. Here is the boundary of A, A(p) = /o(p) is the 
chemical potential associated to the microscopic Hamiltonian 7i, and Aq is the chemical 
potential of the boundary reservoirs. Finally the initial condition for (13. 9p is obtained as 
the limiting empirical density of the chosen microscopic initial configuration of particles. 

We obtain an equilibrium model either if A is the torus and there is no external field 
or in the case of boundary driven systems in which the external field in the bulk matches 
the driving from the boundary; in particular if Aq is constant and E vanishes. In the 
other cases the stationary state supports a non vanishing current and the systems is out 
of equilibrium. 

The coefficients D and x ^i^^e related by the Einstein relation D = R~^Xi where R 
is the compressibility: R^^ = /g , in which /o is the equilibrium free energy associated 
to the Hamiltonian 7i, see [50]. For gradient lattice g defined in Section 12.51 the 

diffusion matrix D and the mobility x are multiples of the identity. For non gradient 



dtp = V- [-D{p)Vp-x{p)E 



(3.9) 



\{p{t,u)) = Xo{u) 



u e dA 



(3.10) 
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models in general D and x ^"^^ iiot diagonal, however, as shown in [51, Lemma 8.3], if the 
Hamiltonian Ti is invariant w.r.t. rotation of 7r/2, then D and x are diagonal. 



We next discuss the large deviation properties of the empirical density; the derivation 
can be found in [3,4,39,50]. Fix a smooth trajectory p = p{t,u), {t,u) G [0,T] x A. We 
want to compute the asymptotic probability that the empirical density tt^ is in a small 
neighborhood of p. If p is not a solution to (13. 9p . this probability will be exponentially 
small and the corresponding rate is called the large deviation dynamical rate functional. 

Consider an initial configuration r] whose empirical measure approximates, as di- 
verges, p(0) and let be the law of the microscopic process starting from such initial 
condition. The dynamical large deviation principle for the empirical density states that 

^ p) ~ exp { - N'lio,T]{p)} (3.11) 

where the rate functional /[o,t] is 

ho,T]ip) = ^£dt{VH,x{p)VH) (3.12) 

in which ( , ) denotes integration in the space variables and VH = VH{t, u) is the extra 
gradient external field needed to produce the fluctuation p, namely such that 

dtp = V- \^-D{p)Vp - x{p){E + VH)] (3.13) 

The interpretation of (13.121) is straightforward; since x is the mobility, /[o,t](p) is the 
work done by the external field Vif to produce the fluctuation p in the time interval 

[o,r]. 

3.2. Thermodynamic functionals and Hamilton- Jacobi equation 

Consider the following physical situation. The system is macroscopically in the sta- 
tionary profile p = p{u), m G A (a stationary solution to (13.91) ) at t = — oo, but at t = 
we find it in the state p. We want to determine the most probable trajectory followed 
in the spontaneous creation of this fluctuation. According to (13.111) this trajectory is the 
one that minimizes / among all trajectories p(t) connecting p to p in the time interval 
[—00,0]. We thus define the so-called quasi-potential as 

Vip) = ^ inf /[_oc,o](p) (3.14) 

p:p(-oo)=p 

m=p 

As shown in [3,5], the functional V solves the Hamilton- Jacobi equation 



1 /-SV , ._6V\ /6V 



2 \ 5p ' 6p / \ 6p 



V^, X(P) V— ) + V ■ -D{p)Vp - x{p)E ) = (3.15) 



2 



note that there is no uniqueness of solutions, e.g. \^ = is always a solution. In [3] it is 
discussed the appropriate selection criterion, that is V is the maximal solution to (I3.15p . 

If the system is in equilibrium then the quasi-potential V coincides with the variation 
of the equilibrium free energy associated to the profile p. The latter can be characterized. 
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by the Einstein fluctuation formula, as the rate of the asymptotic probability of observing 
a given density profile in the equilibrium measure. Namely, if /i at is the invariant measure 
of the generator L^r, then 

/iiv(vr^ ~ p) ~ exp { - N^V{p)} (3.16) 

This relation holds also for nonequilibrium systems, see [3, 16] and, in this sense, the 
solution to the variational problem fl3.14l) is the appropriate generalization of the free 
energy for nonequilibrium systems. Finally, as discussed in [3], for generic nonequilibrium 
models the quasi-potential V is a. non-local functional of p. Notable exceptions are the 
zero-range model and the case, discussed in the Section [3^ of systems with weak external 
field and periodic boundary conditions. We next recall some results on the quasi potential 
for specific lattice gases. 

Zero range process. We consider the zero range process as introduced in Section [2l2] either 
in the torus or in a bounded domain with a weak external field E. Recalling (12.191) and 
(I2.20p . we define the function $ : IR+ IR+ as the activity corresponding to the density 
a, i.e. such that 

a = ———\^k^-i- 3.17 

Z($(a)) ^ g{k)\ ^ ^ 

where Z{ip) is defined in (12.201) . In other words a i— > $(«) is the inverse of the function 
ip I— > R{(p) defined by 

^M = ^f^ (3.18) 

As shown in [2,3, 17,39], the hydrodynamic equation for the zero range process is then 
(13. 9p with D = ^' and x = ^- the case of independent random walks, i.e. g{k) = /c, $ 
is the identity so that D = 1 and x(p) = p. 

Since for the zero range process, as discussed in Section 12.21 the invariant measure 
is always a product measure, the quasi potential ^ is a local functional. Its form can 
be computed directly from the invariant measure by requiring that (I3.16P holds. On the 
other hand it is also possible to solve explicitly the Hamilton- Jacobi equation (13.151) . As 
shown in [2,3] we get 



V{p) = / du 



A 



p(n)log^^^-log^i^M!f))) 
p(nJlog log 



(3.19) 



where (p{u) = $(p(m)) is the stationary activity profile, p being the stationary solution to 
(13. 9p . i.e. the stationary density profile. Equivalently ip solves 

(p{u) = exp{Ao(M)} u e dA ^ ^ 

which, recalling (12. 9p . is just the continuos limit of (I2.18p . 

Boundary driven symmetric exclusion process. We consider here the one-dimensional 
symmetric exclusion process as introduced in Section with A = (0, 1). Let pQ and pi 
be the boundary densities. As shown in [3,26,27], the hydrodynamic equation is (13.90 
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with D = 1 and x(p) = p(l~p)- For this model, if po 7^ Pi the quasi potential is non local, 
which is the signature of macroscopic long range correlations. The quasi potential cannot 
be written in a closed form, but can be obtained by solving a one-dimensional boundary 
value problem. This has been proven in [19,20] by combinatorial methods and in [3,4] by 
the dynamical/variational approach here presented. The result is the following. 

V{p) = snp[ [plog4 + (l-p)logi^ + log— ^1 (3.21) 
/ io L / 1 - / Pi - PoJ 

where the supremum is carried out over all strictly monotone smooth functions / satisfying 
the boundary conditions /(O) = po, /(I) = pi- It has also been shown that there exists a 
unique maximizer for the variational problem (13.211) which is the unique strictly monotone 
solution to the non-linear boundary value problem 

/(I - /)(f)2 + / = P .322) 

/(0)=Po, /(l)=Pi ^ • ' 

in which p = p{u) is the prescribed fluctuation. Knowing that (13.2 ip is the answer, 
the proof amounts to some lengthy but straightforward computations in showing that it 
solves the Hamilton- Jacobi equation (13.151) . see [4] for the details. From (I3.2ip . since V 
is expressed as the supremum of convex functionals we get "for free" that \^ is a convex 
functional. However, as shown below, this convexity property does not hold in general. 

Variational formulae like (13.210 are typical in statistical mechanics, but here the inter- 
pretation it is rather unclear. Firstly it appears strange that we need to maximize and not 
to minimize, secondly the meaning of the test function / is not apparent. For the second 
issue we mention that a dynamical interpretation of / in terms of the hydrodynamics of 
the time reversed process is discussed in [3]. For the first issue we shall show that it is 
connected with the convexity properties of the mobility x- 

We mention that an expression similar to (I3.2ip has also been obtained for the bound- 
ary driven weakly asymmetric exclusion process in [25]. Also the Hamilton- Jacobi ap- 
proach can be applied successfully, see [11]. 

Boundary driven KMP process. We consider here the KMP process introduced in Sec- 
tion 12.41 The hydrodynamic equation is (13.90 with D = 1 and x(p) = P^- Note that 
here p is the energy density and not the particle density as for lattice gases. Similarly to 
the boundary driven symmetric exclusion process, as shown in [12], the quasi potential 
can be obtained by solving a one-dimensional boundary value problem. The result is the 
following. 

\/(p)=inf^(p,/) (3.23) 

where 

^P . , P , /' 







Q{p. /) = / rf« ^ - 1 - log ^ - log jr^^ (3.24) 



and the infimum is carried out over all strictly monotone smooth functions / satisfying 
the boundary conditions /(O) = Tq, /(I) = Ti. It has also been shown that there exists a 
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unique minimizer for the variational problem (I3.23P which is the unique strictly monotone 
solution to the non-linear boundary value problem 

J WY~ J -~P (3.25) 
/(0) = To, /(l)=Ti ^ ' 

in which p = p{u) is the prescribed fluctuation. As the for the boundary driven symmetric 
exclusion process, knowing that (13.231) is the answer, the proof amounts to some lengthy 
but straightforward computations in showing that it solves the Hamilton- Jacobi equation 
(I3.15p . Unlike the boundary driven symmetric exclusion process, the quasi potential for 
the KMP process is not convex. 

A possible interpretation of (I3.23P is the following. The local functional Q{p,f) can 
be thought of as a joint rate functional for both the energy density p and the function /, 
which we can interpreted as a temperature profile. Then the minimization procedure of 
(I3.23P corresponds to the application of a contraction principle. We therefore search for 
the best hidden temperature profile / associated to the energy density profile p. This is 
the inspiring idea behind the following toy model for the invariant measure. 

We will show that the functional V in (I3.23P is the large deviations rate functional of 
a measure on R^^^ which is "simple" enough to be described explicitly and "rich" enough 
to produce such a non-local rate functional. Recall that in Section [2^ we have obtained 
an explicit representation of the invariant measure of the KMP process with a single 
oscillator as a convex combinations of exponential distributions. 

We assume Tq < Ti and let ti, ■ ■ ■ , tAr_i be independent uniform random variables on 
the interval [Tq, Ti]. Denote by < tpj < ■ ■ • < t\^N-i] be order statistics of ti, ■ ■ ■ , t^^i, 
i.e. is the smallest among the tj, t[2] the second smallest and so on. Denote by qn-, 
the distribution of the random vector t[i],- ■■ ,t[Ar_i]; note that qn is a probability on 
[To, Ti]^"-*^. We then define vn as the probability measure on R^'^ whose density w.r.t. 
the Lebesgue measure = YIxgAn given by 



f QNidti,- ■ ■ ,dtN~i) J]^ ^exp{-^a;/ta;} (3.26) 



That is z/tv is a mixture of the exponential Gibbs distribution with temperature profile 
T{x/N) = tx- The measure is not the invariant measure of the KMP process; if we 
compare (13.260 for a single oscillator, N = 2, with the exact expression in (I2.32p we see 
that we replaced the arcsin distribution in [To,Ti] with the uniform one. As diverges, 
the measure un is however a good approximation of the true invariant measure in the 
sense that it leads the rate function in (13.230 . In particular it has the correct asymptotic 
form of the two point correlations. 

To prove the above statement, let us consider the probability measure z/^r on the space 
X [To,Ti]^^ given by 

i>M{d^,dt) = g^idt) Yl ^exp{-^Jt,}dCx (3.27) 

SO that in (I3.26P is obtained as the first marginal of pn, i.e. integrating on the second 
variable t. Recalling the definition (13.60 of the empirical density vr^, we likewise define 
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the empirical temperature profile by requiring that for each smooth function G on 
A {t^ ^G) = ;^ XlxgAjv ^(^/^)^^ Given a smooth function p : A ^ R+ and a smooth 
strictly increasing function / : A ^ [To,Ti] such that /(O) = Tq and /(I) = Ti, we claim 
that 

z>^ (vr^ ^ p , ^ /) ~ exp { - iV^(p, /) } (3.28) 

where Q was defined in (13.241) . To obtain this result, we first observe that if ei, ■ ■ ■ , cat 
are independent exponential random variables with parameter T, then 



1 ^ 



We also recall, see e.g. [28, 1.6], that the random variables Ai := tjij — Tq, A2 : = 
^[2] — ^[1], • • ■ 5 := Ti — t[Af-i] are distributed according the product of A^ exponen- 
tial conditioned on Ai + ■ ■ ■ + Ajv = Ti — Tq. We then get 

QN (r^ ^f)r^exp{-Njyu[- log } 

Since, conditionally on the random variables t^^], x G A^r, the distribution of ^ is the 
product of exponentials, (13.281) follows. Finally, from (13.281) . by maximizing over the 
possible values of /, we easily get that 

^/jvK~p) ~exp{-Arinf6;(p,/)} (3.29) 



3.3. Macroscopic correlation functions 

In Sections 12. 31 and 1 2.41 we found exact formulae for the two-point correlation functions 
Gn{x, y) = Hn{jIx] Vy) of the invariant measure fiN, both for the one dimensional boundary 
driven symmetric exclusion process, see (12.271) . and for the boundary driven KMP process, 
see (I2.38p . For both models we found that, out of equilibrium, they admit long range 
correlations of order 1/A^. More precisely we have. 



CN{x,y) = 




with 

C{u,v) = G{v,u) = -{pi- pofuil-v) , 0<u<v<l (3.30) 
for the boundary driven symmetric exclusion process, and 

G{u,v) = C{v,u) = +{Ti-Tofu{l-v) , 0<u<v<l (3.31) 

for the boundary driven KMP process. Notice that the above functions (13.301) and (13.311) 
only differ by a sign. Moreover the off-diagonal part of the above covariance is proportional 
to the Green function of the Laplacian on the interval [0, 1] with Dirichlet boundary 
conditions, 

A'\u,v) = A^^{v,u) = -u{l-v) , 0<u<v<l (3.32) 
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namely the solution to the problem (9^ A~^(m, v) = 5{u — v), < u, v < 1, with boundary 
condition A~^{u,v) = if either m or f is or 1. 

In this section we will derive the above results from a purely macroscopic point of view. 
More precisely, we consider a one-dimensional boundary driven system with A = (0, 1) and 
no external field and we assume that the transport coefficients in (13. 9p are of the following 
form. The diffusion coefficient is constant, we set D{p) = 1, and x{p) is quadratic so that 
x" is constant. We show that such models have positive, resp. negative, correlations if 
X" > 0, resp. x" < 0- 

Recall that the quasi-potential V{p) solves the Hamilton- Jacobi equation (13.151) . which 
in this context reads 

V^,x(p)V^-Vp)=0 (3.33) 



Sp Sp 

The functional V assumes its minimum at p, the stationary solution to (13. 9p . which in this 
case is a linear function. The correlation function C{u, v), which measures the covariance 
of the density fluctuations with respect to the invariant measure, is then obtained in the 
quadratic approximation of V, i.e. 

V{p) = l{ip- P\ C-\p - p)) + 0((p - pf) (3.34) 

where denotes the inverse operator of C . We can therefore get an equation for C(m, v) 
by expanding the Hamilton- Jacobi equation (13.331) up to second order in {p — p). 

It is convenient to introduce the "pressure" G{h), see [3], defined as the Legendre 
transform of the quasi-potential V{p), 

G{h) = sup \{h,p)-V{p) 

Here h = h{u) can be interpreted as a chemical potential profile. By Legendre dual- 
ity, equation (13.331) can be rewritten as the following Hamilton- Jacobi equation for the 
pressure, 

:^''-^(f)^"-^f)-° f^'^^) 

for any h which satisfies the boundary conditions /i(0) = h{l) = 0. Moreover the expan- 
sion (13.341) gets translated into the following expansion for G, 

G{h) = {h,p) + ^{h,Gh) + 0{h^) (3.36) 

Hence the macroscopic correlation function G can be obtained by expanding equation 
(I3.35P up to second order in h. 
From equation (13.361) we have 

^'^ p{u) + {Gh){u) + 0{h^) (3.37) 



5h{u) 

If we thus plug (I3.37P into (I3.35P and neglect the terms of order we get 

Vh,x{p)Vh-VGh^ = (3.38) 
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for all chemical potential profiles h such that /i(0) = h{l) = 0. To derive the above 
equation we used the fact that p is linear. The macroscopic correlation function C can 
then be determined as the solution to equation f l3.38p satisfying the boundary condition 
C(m, = if m 7^ f and either u or f is or 1. This condition is due to the fact that the 
values of the density at the boundary is fixed by the reservoirs. 

We next define the nonequilibrium contribution to the covariance as the function B 
such that 

C{u,v) = x{.p{.u))5{u-v) + B{u,v) , u,veA (3.39) 

Note that, since D = 1, is the local equilibrium variance. By plugging (13.391) into 

(I3.38p . we get that B solves 

{dl + dl)B{u, v) = -{Vpf x" S{u - v) (3.40) 

together with the boundary condition B{u,v) = if either m or f is or 1. The above 
equation can also be derived within the fluctuating hydrodynamic theory, see [49]. Hence 

B{u,v) = -^iVpfx"A~\u,v) 

which, by (13.321) and recalling that x{p) = pi)-— P) for the exclusion process and x(p) = P^ 
for the KMP process, agrees with (I3.30p and (13.311) . 

In [10] we derive the equation satisfied by the off diagonal covariance B for arbitrary 
dimension, D, Xi ci^nd external field E. This equation allows to establish, for a class of 
models, whether the correlations are positive or negative. 

3.4- Weakly asymmetric models with periodic boundary conditions 

We consider here a lattice gas with periodic boundary conditions, namely A is the 
(i-dimensional torus, and constant weak external field E. As discussed in Section 12. 5[ 
from a microscopic point of view, if the model is gradient then the invariant measure does 
not depend on the external field E. As we show here, from a macroscopic point of view, 
any system behaves as gradient models. 

The precise statement is the following. Consider the variational problem (I3.14p defin- 
ing the quasi-potential V in the present setting of periodic boundary conditions and 
constant external field E. Then V does not depend on E and therefore coincides with the 
solution to (I3.14P with = 0, namely with the free energy associated to the microscopic 
Hamiltonian Ti. 

We suppose given the transport coefficients D and x iii (I3.12p - (l3.13p so that the 
Einstein relationship D{p) = R[p)~^x{p) holds; recall that while D and x ^^e matrices, 
the compressibility i? is a scalar. In the case of periodic boundary conditions and constant 
field E there is a one parameter family of stationary solutions to (13.90 which are simply 
the constant functions p{u) = m, m & M+. Given m G M+ we define 

which is a strictly convex function with minimum at p = m. We claim that the solution 
of the variational problem (I3.14p with p = m is the functional 

J'M = [ dufMu)) (3.41) 

J A 
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for any value of the external field E. 

li E = 0, by using the Einstein relation D{p) = fm{p)x{p)^ it is easy to check that JF„ 
solves the Hamilton- Jacobi (13.151) . If is a constant, since the boundary conditions are 
periodic, we have that 

<^,V.x(p)i^>=0 

hence JF^ solves the Hamilton-Jacobi f l3.15p for any (constant) external field E. It is also 
not difficult to check that JF„^ is the maximal solution to the Hamilton-Jacobi equation 
(I3.15P so that the claim is proven. 

3.5. Glauber + Kawasaki 

We consider here the macroscopic behavior of the Glauber -|- Kawasaki process intro- 
duced in Section F2.6[ The empirical density is defined as in (13.61) . We emphasize that 
in this model the empirical density is not locally conserved due the reaction terms in 
the microscopic dynamic (I2.48p . Accordingly, the hydrodynamic equation is given by the 
reaction diffusion equation 

dtp = ^Ap + bip) - dip) (3.42) 

where the reaction terms b and d, which are polynomials in p, are determined by the rates 
Cxiv) in fl^^iH]) as follows, [18,35], 

Kp) = i^p(co(^)(l - ^o)) , d{p) = iyp{co{r])rio) (3.43) 

where Up is the Bernoulli measure with density p. In particular, in the reversible case 
where the rates Cx{ri) are as in (12.520 . b{p) and d{p) have the form 

b{p)=A,{l-p)^{p) , d{p) = A2pM (3.44) 

where (p{p) is the expected value of h{T]) in (12.520 with respect to Up and Ai, A2 > 0. 
We consider this system only with periodic boundary conditions. The equilibrium profile 
thus corresponds to a constant density p which solves b{p) = d{p) and gives an absolute 
minimum of the potential U, defined by U'{p) = —[b{p) — d{p)]. 

The associated large deviation asymptotics is in the same form as in (13. lip , but here 
the rate functional /[o,t] is not simply quadratic in the external field. Indeed in [35] it is 
proven that it is given by 

I[o,T]{p) = £dt[^{VH,p{l~p)VH) 

+ (b{p), (1 - + He'') ) + (d(p), (1 - e-^ - i^e"^))} (3.45) 

where the external potential H is connected to the fluctuation p by 

dtp = ^Ap - V ■ (p(l - p)VH) + b{p)e'' - d{p)e-'' (3.46) 

As in Section ing we analyze the variational problem (I3.14p . The associated Hamilton- 
Jacobi equation [1] is 

fl(p,— )=0 (3.47) 
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where the "Hamiltonian" Sj is not anymore quadratic in the momenta and it is given by 

i3(p, H) = ^{H, Ap) + ^{VH, p(l - p)VH) - (6(p), 1 - e^> - (dip), 1 - e"^) (3.48) 

If b and d are as in (13.441) it is easy to find the solution V of (13.471) . [32]. Let p = 
Ai/ {Ai + A2), the unique root of b{p) — d{p) = 0, then 

V{p) = / du p log ^ + (1 - p) log —-^ (3.49) 
io L p 1-pJ 

If the reaction rates Cx{ri) are of the form (I2.52p . then the invariant measure is Bernoulli 
and (I3.49P follows. On the other hand, as shown in [32], there are choices of the reaction 
rates such that (12.521) fails but (I3.44p holds. In this cases (I3.49P still holds and we may 
say that reversibility is restored at the macroscopic level or that time reversal invariance 
is violated "weakly" by the microscopic dynamics. 

Correlation functions. In Section we studied long range correlations for some bound- 
ary driven (hence non equilibrium) conservative models. Here we consider equilibrium 
states for the Glauber + Kawasaki dynamics, which is non conservative, and we study 
their macroscopic correlation functions. In particular we show that, if the microscopic 
dynamics violates time reversal invariance "strongly", that is (I3.44p does not hold, long 
range correlations do appear, [1]. 

Recall that, in order for the system to be reversible, the rates Cx{f]) of the Glauber 
dynamics should be of the form (12.521) . Their relationship with the coefficients b{p) and 
dip) in ([332]) is given in (K43\\ . 

Equation (I3.47P is a very complicated functional derivative equation which, as in Sec- 
tion \'S.'S\ can be solved by successive approximations by formal power series expansion in 
p — p. Here p is a constant stationary solution of (I3.42p . i.e. a root of 6(p) — d{p) = 0. 

Let C be the correlation function. It gives the second order approximation of the 
quasi-potential as in (13.341) . Proceeding as in Section [3^31 see [1] for further details, we 
get that C solves 

^dlC{u, v) - {d, - b,)C{u, v)-^ p(l - p) dl5{u -v) + bo6{u - v) = (3.50) 

where 

bi = b'{p) , di = d'{p) , 60 = Kp) = d{p) = do 
Notice that, if (imj) holds, we get 

7:=feo-p(l-p)(rfi-6i) = 0. (3.51) 

and in this case, recalling (I3.49p . we of course have C{u, v) = p(l — p) 6{u — v). Conversely, 
a solution of the form C{u, v) = a6{u — v), for some a > 0, exists only if (13.511) holds and 
therefore a = p(l — p). 

The above considerations imply that long range correlations do appear whenever (13.511) 
fails. In this case we say that irreversibility persists at the macroscopic level. As in Section 
13.31 we introduce the off-diagonal covariance B such that 

C{u, v) = p{l - p)5{u - v) + B{u, v) 
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and we get that B solves 



- ^dlB{u, v) + {di - bi)B{u, v) = -fS{u - v) (3.52) 

where 7 is defined in (13.511) . Note that di — bi, being the second derivative of the potential 
calculated in a minimum, is positive. Let R = [di — bi — (1/2) A)~^ be the resolvent of 
the Laplacian on the torus. Then the solution of (13.521) is 

B{u,v) =-f R{u,v) (3.53) 

Since R{u,v) > 0, we conclude that the correlation B{u,v) has the same sign as 7. 

While 7 = corresponds to the macroscopically reversible situation, in general 7 may 
have either sign. For instance, given a G (—1,1), take the fiip rates given by 

coiv) = vo[^- " — 2 — j + ~ + " — 2 — / 

for a > presence of surrounding particles enhances the birth rate and suppresses the 
death rate. We thus expect that the two point correlation to be positive for a > and 
negative for a < 0. We have 7 = a{l — a){2 — a)~^ which shows that this is indeed the 
case. 

In [18] it is shown that fiuctuations from the hydrodynamical equation with standard 
Gaussian normalization converge, as ^ 00, to an Ornstein-Uhlenbeck process. The 
stationary correlations of this process agree, as they should, with (13.531) . 

3. 6. Boundary driven asymmetric exclusion process 

We start from the representation of the invariant measure for the boundary driven 
totally asymmetric exclusion process obtained in [23] and illustrated in Section 12.71 We 
call vr^ the empirical measure associated to the configuration 77, 7^ the empirical measure 
associated to ^ and ^(p, /) the joint rate functional 

z/7v(7r^^P,7^^/) -e-^^^''-^) (3.54) 

where un is the measure (12.561) . From formula (I2.57p . using the contraction principle of 
large deviations, we obtain directly 

Piv (vr^ ^ p) ~ e-^^^") (3.55) 

where 

y(p)=inf^(p,/) (3.56) 

This argument suggests a different representation, from the one obtained in [21], for the 
non local rate functional V{p). In particular, while in [21] V{p) is obtained either as an 
infimum or a supremum, depending on the values of the chemical potentials Aq and Ai, 
here we write V{p) always as an infimum. 

We construct explicitly this new representation in the special case Aq = Ai = 0. In this 
case the measure is uniform on the set of complete configurations, defined by (I2.55p . 
and the joint rate functional Q is easily obtained as a restriction of the one associated to 
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the uniform measure over all configurations (?7,0 ^ -^^^ ^ X^^ . We define the set of 
complete profiles 

(p, /) : fdv [p{v) + fiv)] >u, ue [0, 1]; f'dv [p{v) + f{v)] = l] 



C :-- 

Remember that p and / are density profiles for configurations of particles satisfying an 
exclusion rule so that they take values in [0, 1]. Then in the special case Aq = Ai = we 
have 



fdu[h{p) + h{f)\ if (p,/)eC 



<3{pJ)={ rv/^^ - - vH.j,^^ (3.57) 

-oo otherwise 

where h{x) = xlog(2x) + (1 — x) log [2(1 — x)]. Note that we do not need to add a 
normalization constant due to the fact that the constant profiles (|, |) belong to C and 
Q (f ) ^) = 0. Using (13.561) we obtain the following variational representation for the 
quasi-potential V, 

V(p) = inf [ du \h{p) + hif)] (3.58) 
that has to be compared with the one in [21] 

V{p)= sup / dM{plog[p(l-/)] + (l-p)log[(l-p)/]+log4} (3.59) 

{/:/{0)=l,/(l)=0}7o 

where the supremum is over monotone functions. 

In [21] it is shown that the supremum in (I3.59P is obtained when f = F'^ and 



F,{u)=CE(^j^ dv[l-p{v)] 



where CE means Concave Envelope. We will show next that F'^ is also the minimizer of 
the problem (I3.58p . This is equivalent to prove that 

inf ! duh(F'(u))= [ duhiF'iu)) (3.60) 

Note that F' has to be a density profile so that F is increasing. Moreover F is defined 
up to an additive constant so we can choose -F(O) = 0. The condition (p, F') G C easily 
reads as 



F(0) = 0, F{u)> / dv[l- p{v)]yue [0,1] , F{1)= / dv[l 
Jo Jo 



p{v)] (3.61) 



It is clear that if F satisfies condition (13.611) . then also CE{F) satisfies (I3.6ip : or, equiv- 
alently, if (p, F') G C then also (p, CE{F)') G C. Moreover the following elementary 
inequality holds due to the convexity of /i, 

fdu h{F\u)) >{h- a)h(^^^-^^^) (3.62) 
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which immediately imphes 



duh{F'{u)) > / duh{[CE{F)]'{u)) 
Jo 

We thus conclude that we can restrict the infimum fl3.60p over the set of concave func- 
tions F satisfying conditions fl3.6ip . Still a direct application of (13.621) imposes that the 
minimizer has to be the smallest among them, that is Fp. 

Using the above result we can finally prove the equivalence between the two different 
representations (13.581) and (I3.59p of V{p). In order to prove it we just have to show that, 
for any density profile p, 

f du {p\og{l-F'^) + {l-p) log F'^}= rdn{(l-F;)log(l-F;) + F;iogF;} (3.63) 
Jo Jo 

The contributions on both sides in (I3.63P from the domain of integration where F^ = 1 — p 
are clearly equal. Consider now a maximal interval [a, b] where F'^ ^ 1 — p. Then on this 
interval we have ^ 

From this fact, an easy computation shows that also the contributions on both sides of 
from the integrations over [a, h] are equal. 



Acknowledgments 

We are grateful to A. Faggionato, G. Basile and J. Lebowitz for very illuminating discussions 
and collaborations. G. J.-L. would like to thank the organizers of the program Principles of the 
dynamics of nan equilibrium systems held at the Newton Institute in 2006, for the invitation 
to participate. Thanks in particular to C. Godreche for suggesting to write this paper. We 
acknowledge the support of PRIN MIUR 2004-028108 and 2004-015228. 



References 

[1] Basile G., Jona-Lasinio G., Equilibrium states with macroscopic correlations. Inter. J. Mod. 
Phys. B 18, 479-485 (2004). 

[2] Bertini L., De Sole A., Gabrielli D., Jona-Lasinio G., Landim C., Fluctuations in stationary 
non equilibrium states of irreversible processes. Phys. Rev. Lett. 87, 040601 (2001). 

[3] Bertini L., De Sole A., Gabrielli D., Jona-Lasinio G., Landim C., Macroscopic fluctuation 
theory for stationary non equilibrium state. J. Statist. Phys. 107, 635-675 (2002). 

[4] Bertini L., De Sole A., Gabrielli D., Jona-Lasinio G., Landim C., Large deviations for the 
boundary driven simple exclusion process. Math. Phys. Anal. Geom. 6, 231-267 (2003). 

[5] Bertini L., De Sole A., Gabrielli D., Jona-Lasinio G., Landim C., Minimum dissipation 
principle in stationary non equilibrium states. J. Statist. Phys. 116, 831-841. (2004) 



32 



[6] Bertini L., De Sole A., Gabrielli D., Jona-Lasinio G., Landim C, Current fluctuations in 
stochastic lattice gases. Phys. Rev. Lett. 94, 030601 (2005). 

[7] Bertini L., De Sole A., Gabrielli D., Jona-Lasinio G., Landim C, Non equilibrium current 
fluctuations in stochastic lattice gases. J. Statist. Phys. 123 237-276 (2006). 

[8] Bertini L., De Sole A., Gabrielli D., Jona-Lasinio G., Landim C, Large deviation approach 
to non equilibrium processes in stochastic lattice gases. Bull. Braz. Math. Soc, New Series 
37 611-643 (2006). 

[9] Bertini L., De Sole A., Gabrielli D., Jona-Lasinio G., Landim C, Large deviations of the 
empirical current in interacting particle systems. Theory of Probability and its Applica- 
tions, 51, 2-27 (2007). 

[10] Bertini L., De Sole A., Gabrielli D., Jona-Lasinio G., Landim C., On the long range corre- 
lations of thermodynamic systems out of equilibrium. In preparation. 

[11] Bertini L., Gabrielli D., Landim C., Large deviations for the weakly asymmetric boundary 
driven exclusion process. In preparation. 

[12] Bertini L., Gabrielli D., Lebowitz J.L., Large deviations for a stochastic model of heat flow. 
J. Statist. Phys. 121, 843-885 (2005). 

[13] Bodineau T., Derrida B., Current fluctuations in non-equilibrium diffusive systems: an 
additivity principle. Phys. Rev. Lett. 92, 180601 (2004). 

[14] Bodineau T., Derrida B., Distribution of current in non-equilibrium diffusive systems and 
phase transitions. Phys. Rev. E 72, 066110 (2005). 

[15] Bodineau T., Derrida B., Current large deviations for asymmetric exclusion processes with 
open boundaries. J. Stat. Phys. 123, 277-300 (2006). 

[16] Bodineau T., Giacomin G., From dynamic to static large deviations in boundary driven 
exclusion particles systems. Stoch. Proc. Appl. 110, 67-81 (2004). 

[17] De Masi A., Ferrari P., A remark on the hydrodynamics of the zero-range processes. J. 
Statist. Phys. 36, 81-87 (1984). 

[18] De Masi A., Ferrari P., Lebowitz J.L., Reaction- Diffusion Equations for Interacting Particle 
Systems. J. Statist. Phys., 44, 589-644, (1986). 

[19] Derrida B., Lebowitz J.L., Speer E.R., Free energy functional for nonequilibrium systems: 
an exactly solvable case. Phys. Rev. Lett. 87, 150601 (2001). 

[20] Derrida B., Lebowitz J.L., Speer E.R., Large deviation of the density profile in the steady 
state of the open symmetric simple exclusion process. J. Statist. Phys. 107, 599-634 (2002). 

[21] Derrida B., Lebowitz J.L., Speer E.R., Exact large deviation functional of a stationary open 
driven diffusive system: the asymmetric exclusion process. J. Statist. Phys. 110, 775-810 
(2003). 

[22] Dorfman J.R., Kirkpatrick T.R., Sengers J.V., Generic long-range correlations in molecular 
fluids. Annu. Rev. Phys. Chem. 45, 213-239 (1994). 



33 



[23] Duchi E., SchaefFer G., A combinatorial approach to jumping particles. J. Combin. Theory 
Ser. A 110, 1-29 (2005). 

[24] Einstein A., Theorie der Opaleszenz von homogenen Fliissigkeiten und Fliissigkeitsgemis- 
chen in der Ndhe des kritischen Zustandes. Annalen der Physik, 33, 1275-1298 (1910), 
English translation in The collected papers of Albert Einstein, vol.3 p. 231-249, Princeton 
University Press, 1993. 

[25] Enaud C, Derrida B., Large deviation functional of the weakly asymmetric exclusion pro- 
cess. J. Stat. Phys. 114, 537-562 (2004). 

[26] Eyink G., Lebowitz J.L., Spohn H., Hydrodynamics of stationary nonequilibrium states for 
some lattice gas models. Commun. Math. Phys. 132, 253-283 (1990). 

[27] Eyink G., Lebowitz J.L., Spohn H., Lattice gas models in contact with stochastic reservoirs: 
local equilibrium and relaxation to the steady state. Commun. Math. Phys. 140, 119-131 
(1991). 

[28] Feller W., An introduction to probability theory and its applications, Volume 2. Wiley, 1966. 

[29] Freidlin M.I., Wentzell A.D., Random perturbations of dynamical systems. Springer, 1998. 

[30] D. Gabrielli, G. Jona-Lasinio, C. Landim, Onsager reciprocity relations without micro- 
scopic reversibility. Phys. Rev. Lett. 77 (1996), 1202-1205. 

[31] D. Gabrielli, G. Jona-Lasinio, C. Landim, Onsager symmetry from microscopic TP invari- 
ance. J. Statist. Phys. 96 (1999), 639-652. 

[32] Gabrielli D., Jona-Lasinio G., Landim C., Vares M. E., Microscopic reversibility and ther- 
modynamic fluctuations. In "Boltzmann's Legacy 150 Years after his Birth" , Atti dei Con- 
vegni Lincei, 131, 79-87 (1997). 

[33] Giardina C., Kurchan J., Redig F., Duality and exact correlations for a model of heat 
conduction. Preprint 2006, |cond-mat /06 1 2 198 1 

[34] Harris R.J., Schiitz G.M., Fluctuation theorems for stochastic dynamics. Preprint 2007, 
|cond-mat/0702553 . 

[35] Jona-Lasinio G., Landim C., Vares M.E., Large deviations for a reaction-diffusion model. 
Probab. Theory Related Fields 97, 339-361 (1993). 

[36] Katz S., Lebowitz J.L., Spohn H., Nonequilibrium steady states of stochastic lattice gas 
models of fast ionic conductors. J. Statist. Phys. 34, 497-537 (1984). 

[37] Kingman J. E.G., Markov population processes. J. Appl. Prob. 6, 1-18 (1969). 

[38] Kipnis C., Marchioro G., Presutti E., Heat flow in an exactly solvable model. J. Statist. 
Phys. 27, 65-74 (1982). 

[39] Kipnis C., Landim G., Scaling limits of interacting particle systems. Springer-Verlag, Berlin, 
1999. 

[40] Landau L., Lifshitz E., Physique Statistique, MIR Moscou 1967. 



34 



[41] Lanford O.E., Entropy and equilibrium states in classical statistical mechanics. Lecture 
Notes in Physics 20, Springer, Berlin, 1973. 

[42] Liggett T.M., Stochastic Interacting Systems: Contact, Voter and Exclusion Processes. 
Grundlehren der mathematischen Wissenschaften 324, Springer, 1999. 

[43] Lovasz L., Discrete analytic functions: an exposition. Surveys in differential geometry. Vol. 
IX, 241-273, Surv. Differ. Geom., IX, Int. Press, Somerville, MA, 2004. 

[44] Maes C, Netocny K., Shergelashvili B., A selection of nonequilibrium issues. Lecture 
notes from the 5th Prague Summer School on Mathematical Statistical Mechanics (2006), 
math-ph/0701047[ 

[45] Onsager L., Reciprocal relations in irreversible processes. I. Phys. Rev. 37, 405-426 (1931). 
//. Phys. Rev. 38, 2265-2279 (1931). 

[46] Onsager L., Machlup S., Fluctuations and irreversible processes. Phys. Rev. 91 (1953), 
1505-1512; Phys. Rev. 91 (1953), 1512-1515. 

[47] Santos J.E., Schiitz G.M., Exact time- dependent correlation functions for the symmetric 
exclusion process with open boundaries. Phys. Rev. E 64, 036107 (2001). 

[48] Schiitz G.M., Exactly solvable models for many-body systems far from equilibrium. In "Phase 
transitions and critical phenomena" vol. 19, C. Domb and J. Lebowitz eds., Academic Press, 
2001. 

[49] Spohn H., Long range correlations for stochastic lattice gases in a nonequilihrium steady 
state. J. Phys. A 16, 4275-4291 (1983). 

[50] Spohn H., Large scale dynamics of interacting particles. Berlin, Springer- Verlag, 1991. 

[51] Varadhan S.R.S., Yau H.T., Diffusive limit of lattice gas with mixing conditions. Asian J. 
Math. 1, no. 4, 623-678 (1997). 



35 



